Optical bistability in subwavelength apertures containing nonlinear media 
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We develop a self-consistent method to study the optical response of metallic gratings with non- 
linear media embedded within their subwavelength slits. An optical Kerr nonlinearity is considered. 
Due to the large E-fields associated with the excitation of the transmission resonances appearing 
in this type of structures, moderate incoming fluxes result in drastic changes in the transmission 
spectra. Importantly, optical bistability is obtained for certain ranges of both flux and wavelength. 
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Since the appearance of the photonic crystal concept, 
there has been a strong interest in the optical properties 
of nano- and micro- structured systems. This is due, in 
part, to the potential applications for small all-optical 
devices. An interesting possibility is to include some 
non-linear elements in the structures^ which may re- 
sult in optical switches and gates. Systems based on 
nonlinear photonic crystals presenting optical bistabil- 
ity have been recently proposed^i^*Sii. Another promis- 
ing route for nonlinear optics is to use structured metal 
films. The enhanced optical transmission phenomenon 
found in subwavelength hole arrays^ and the high local 
field enhancements associated with it-S., suggest the possi- 
bility of strong non- linear transmission effects if nonlinear 
media are embedded in the structure. Recently, photon 
transmission gated by light of a different wavelength were 
reported^*', and the possibility of bistability in enhanced 
transmission has been mentionedii. However, up to our 
knowledge, a detailed calculation for such metallic struc- 
tures containing nonlinear media had not yet been done. 

In this paper, we present a theoretical analysis of 
the transmission properties of a one-dimensional metallic 
grating with subwavelength slit3i2iiiiiiiS4i&, filled with a 
Kerr nonlinear media. As we will discuss below, even this 
simple system presents interesting optical properties, as 
optical bistability for certain ranges of wavelength and 
incident flux. 

In Fig. n we show a schematic view of the structure 
under study. The chosen parameters of the grating are: 
d = 0.75 /im, a = 0.05 //m, and h = 0.45 ^m. This set 
of parameters is typical for experimental studies of sub- 
wavelength apertures in the optical regime^^.. The slits 
are supposed to be filled with a Kerr nonlinear media, 
whose dielectric constant at point r depends on the in- 
tensity of the E-field at this point, |i?(r)p. 



(1) 



where e; is the value of the dielectric constant at low 
intensity, and x'"^"' is the third order susceptibility related 
to the Kerr effect. Since the change of the dielectric con- 
stant is not very large, it is common to approximate the 




FIG. 1: Schematic view of the metallic structure under study: 
a metallic grating of period d, film thickness h, and slit width 
a. Inside the slits, we consider a Kerr nonlinear material, 
characterized by a Kerr coefflcient n2- Light is considered to 
impinge at normal incidence with an E-field pointing across 
the slits. 



refraction index as: n{I) = n{0) + 712/, where / is the 
intensity of the optical field (measured in units of energy 
flux), and 712 is known as the Kerr coefficient. This allows 
to express the incident fiux, Iq, in units of n^^ and the 
results will be valid for all Kerr materials with the same 
n(0). In what follows we consider n(0) = 1.5 (e;=2.25). 

The electromagnetic (EM) properties of these grat- 
ings are analyzed by a modal expansion of the electric 
and magnetic fields. Two simplifications are incorpo- 
rated to the exact modal expansion: we consider per- 
fect metal boundary conditions and, as the slit width 
is much smaller than the wavelength, only the funda- 
mental eigenmode in the modal expansion of the EM 
fields inside the slits is included. We have demonstrated 
in previous theoretical worksiS that these are reason- 
able approximations when analyzing optical properties of 
nanostructured good metals like silver and gold. More- 
over, within this framework, we have been able to re- 
produce in semi-quantitative terms the phenomena of 
beaming^^ and enhanced transmission^*^, appearing for 
subwavelength apertures in the optical regime. We con- 
sider p-polarized light (E-field pointing across the slits), 
for which both extraordinary transmission and electric 
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field enhancements have been reported. Furthermore, in 
tliis paper we concentrate on light impinging at normal 
incidence onto the structure. 

In order to account for the nonlinear response of 
the material within the slits, an iterative self-consistent 
method is used. For this, the slit region is divided into 
TV narrow slices, perpendicular to the z-direction. For 
a given iteration, the magnetic field in the n-th slice is 
expressed as 



Hy - 



Bne 



(2) 



where fcp is the wavenumber of the incident light (in vac- 
cuum), e„ is the dielectric constant in the n-th slice, 
and An and B„ are the modal expansion coefficients. 
Within each slice, the electric field can be obtained from 
Ex = {—i/uje€Q){dHy/dz). By matching Ex and Hy, the 
coefficients of the modal expansion in two consecutive 
strips can be related as follows, 
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with the definitions C± = 1 ± e„+i/e„, '^n,± — 
exp [i(y^± ^e„+i)fcoz„] , and z„ marking the position 
of the interface between slices n and n -|- 1. The modal 
coefficients at the entrance and exit of the slit are re- 
lated through the product of the matrices corresponding 
to the matching of EM fields between consecutive slices. 
Finally, the connection with the incoming and outgo- 
ing waves outside the grating is established, obtaining 
the transmission and reflection coefficients for the whole 
structure. In the next step of the iterative procedure, a 
new set of is obtained through Eq.(l) from the previ- 
ous calculation of the E-fields in the slices. For a given 
value of N and /q, this linear response calculation is re- 
peated until convergency in the set of e„ is reached. And 
then, for each value of /q, N is increased until the trans- 
mitted flux does not depend on N . Still, the described 
procedure does not completely specify the solution: this 
non-linear system may present a multiplicity of solutions, 
and the transmitted flux obtained at the end may depend 
on the profile of e{z) chosen as seed of the iterative pro- 
cedure. Two choices are used in this paper: for a given 
incident flux we use as seeds of the iterative procedure 
the converged dielectric constant profiles obtained from a 
previous calculation for either smaller or larger incident 
fluxes. 

Fig- EI a-) renders the transmittance (total transmitted 
flux divided by incident flux) versus wavelength, for dif- 
ferent values of /q in units of . For the chosen param- 
eters, in the linear regime (red curves in both panels (a) 
and (b)), the metallic grating presents two transmission 
resonances, appearing at A = 1.46 /um and A — 0.78 /im, 
very close to the period of the grating. These two peaks 
are representative of the two types of transmission reso- 
nances appearing in metallic gratingsi^ti^: i) slit waveg- 
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FIG. 2: Transmittance versus wavelength for different values 
of incident flux obtained while increasing (a) and decreasing 
(b) the incident flux. The parameters of the metallic gratings 
are: d — 0.75/im, a — 0.05/xm, and h = 0.45/xm. The incident 
fluxes considered are (1) 10"*n^\ (2) 10~'^n~^ , (3) 10"^nj\ 
and (4) IQ-^n^^ 



uide modes in which the EM-ficlds are highly concen- 
trated within the slits and ii) coupled surface plasmon 
polaritons (SPP) excited at the entrance and exit inter- 
faces of the structure. In Fig. |2Ia), the transmittance 
at a given wavelength has been obtained by consider- 
ing the converged situation at a lower incident flux (for 
that particular wavelength) as the seed for the iterative 
scheme. As can be seen, as the incident flux is increased, 
the peaks shift to larger wavelengths and their magni- 
tude decreases. In Fig.EJb) we represent the same quan- 
tity, transmittance versus wavelength, for the same struc- 
ture, but for decreasing incident fluxes. For a given wave- 



length, the flux is first increased up to /„ 
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which is a representative value of a "high-flux" situa- 
tion, and then decreased slowly to the value /q, using as 
seed for the iterative procedure the converged solution 
obtained for a higher flux. Remarkably, the transmission 
spectra differ for increasing and decreasing fluxes, which 
is a clear manifestation of optical bistability. Besides, 
for decreasing incident flux, the transmission spectra de- 
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FIG. 3: Transmitted flux versus incident flux for wavelength 
(a) A = 1.55/im and (b) A = 0.8/im for a metaUic grating with 
parameters: d — 0.75/im, a = 0.05/im, and h = 0.45/im. Solid 
(dashed) lines correspond to increasing (decreasing) fluxes. In 
the inset of panel (a) it is shown the transmitted flux versus 
increasing incident flux for A = l.bbfim for a metallic grating 
of the same a and d than in the panels but with h = 13/im. 



pends on the initial incident flux /max- In particular, 
in Fig. |2Jb) the transmittance peaks for incident flux 
10~^n^^ would increase (and even reach unity) for larger 

values of /max- 

Optical bistability is better analyzed considering the 
situation for a fixed wavelength. In Fig. |3|[a), we show 
the transmitted flux as a function of Iq for A = 1.55/ito, 
which is the typical wavelength used in telecommuni- 
cations (and slightly larger than the wavelength of the 
waveguide mode in the structure under consideration in 
this paper). A clear bistable loop is observed at a range 
of incident fluxes. When increasing the incident flux (full 
line in Fig. |3Ia)), the transmitted flux increases. At 
a certain incident flux (for this particular case around 
5.4 X 10^^/712), the transmitted flux jumps to a higher 
value in a discontinuous manner. By contrast, the sys- 
tem has a different behavior when the incident flux is 
decreased. The system, now in a high transmittance sit- 
uation, instead of jumping back to low transmitted flux, 
follows the upper branch in Fig. |2Ia). In this system, 
where losses are not included, the upper branch finishes 
when the ratio between the transmitted flux and the in- 
cident flux reaches the unity. Then, the transmitted flux 
jumps discontinuously to a lower value. The presence of 
a narrow bistability region such as the one illustrated in 



Fig. Ola) is an interesting situation for device applica- 
tions, as small changes in the incident flux would imply 
significant changes in the transmitted flux. 

Optical bistability is also obtained associated with the 
transmission resonance linked to the excitation of two 
coupled SPPs. In Fig. |2Ib), we show the bistable loop 
obtained for a wavelength of A = 0.8/im, just slightly 
larger than the period of the grating. The bistable region 
is wider than in the previous case, mainly due to the 
higher incident flux needed to jump to the upper branch 
when increasing the incident flux. 

It is interesting to note here the similarities between 
the curves rendered in Fig. 3 for increasing Iq and the 
measurements of light tunneling through individual pin- 
holes in a granular gold film covered with a layer of a 
nonlinear material (polydiacetylene) reported in Ref. 21. 
These results were explained in terms of a "photon block- 
ade" phenomenon. In our system, if the thickness of the 
metal film is large enough, the spectral separation be- 
tween different slit waveguide modes is strongly reduced. 
Then, as the incident flux is increased, several resonances 
can participate in the transmission process resulting in 
a staircase behaviour of the transmitted flux versus the 
incident one, as can be seen in the inset of Fig. 3a for 
a metallic grating oi h = 13/im. This behaviour is a 
clear signal of the existence of "photon blockade" phe- 
nomena in the system analyzed here. We believe that 
the structure proposed in this paper (array of subwave- 
length apertures) is an alternative that could allow the 
investigation of this interesting phenomenon in a more 
controlled way. 

We find that there are two characteristic EM-field pat- 
terns associated, respectively, with the lower and the 
upper branch of the bistability loop. At a given wave- 
length, the field pattern within the slits for the lower 
branch corresponds approximately to the field pattern 
obtained in the linear response calculation. Even in this 
case, the lower branch presents, close to its instability, 
non-linearities. However, these non-linearities are asso- 
ciated more to the change in the average dielectric con- 
stant inside the slit (via Eq.(^)) than to changes in the 
EM-ficld profile. The upper branch presents a different 
field pattern which is, essentially, the one associated to 
the closest transmission resonance located at a shorter 
wavelengths^. 

Let us now study how the width of the bistability re- 
gion depends on the wavelength. In Fig. 0] we show the 
incident flux corresponding to the start and the end of 
this region as a function of the wavelength. For wave- 
lengths in the interval between 0.79 ^va and 1.53 /xm, 
bistability is related to the shift to larger wavelengths 
of the coupled SPP resonance. The width of the bista- 
bility region increases dramatically with the wavelength. 
We arbitrarily terminate the curves at A=1.3/im, since 
the incident flux needed to jump to the upper branch, 
which corresponds to the end of the bistability region, 
becomes so high that our model is not valid any longer. 
For wavelengths larger than 1.53 /im, the first bistability 
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FIG. 4: Incident flux, in logarithmic scale, corresponding to 
the beginning (continuous line) and end (dashed line) of the 

first bistable region as a function of the wavelength for a 
metallic grating with parameters: d = 0.75/im, a = 0.05/im, 
and h = 0.45/im. 

region is related to the shift to larger wavelengths of the 
slit waveguide mode appearing approximately at 1.46 fim 
at low incident flitx. Notice that the bistability regions 
do not start exactly at the position of the linear-response 
transmission resonances. As an approximate phenomeno- 
logical rule we have found that bistability regions start 
at slightly larger wavelengths, for which transmittance is 



of the order of 1 /4 of the transmission peak. 

Given available fluxes and Kerr nonlinear media, op- 
tical bistability should be experimentally accessible al- 
ready. For instance, some polymers can present a Kerr 
coefRcient up to 10~^ — 10~^cm'^/W. These values would 
imply incident fluxes of the order of 10^ — 10^ W/cm^ for 
observing bistability, for the parameters chosen in this 
article. For higher ratios between the lattice parameter 
and the slit, given that the energy flux inside the slits in- 
creases as d/a, bistability should occur for even smaller 
values of the incident flux. 

In summary, we have analyzed the optical response of a 
nicitallic grating with subwavelength slits containing Kerr 
nonlinear media by means of a self-consistent method. 
The nonlinear response induces changes in the transmis- 
sion spectra, with shifts of the transmission peaks to 
larger wavelengths. In addition, transmission spectra dif- 
fer for increasing and decreasing intensity. An important 
result is the observation of optical bistability for certain 
ranges of wavelength and incident flux. The levels of op- 
tical power needed for observing these effects are achiev- 
able for typical Kerr nonlinear materials. 
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